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The energy spectrum of neutrino-induced upward-going muons in MACRO was analysed in terms of relativity 
principles violating effects, keeping standard mass-induced atmospheric neutrino oscillations as the dominant 
source of — > v T transitions. The data disfavor these possibilities even at a sub-dominant level; stringent 
90% C.L. limits are placed on the Lorentz invariance violation parameter |Au| < 6 x 10~ 24 at sin 26 v = 
and | At; | < 2.5 -j- 5 x 1CP 26 at sin 28 v = ±1. The limits can be re-interpreted as bounds on the Equivalence 
Principle violation parameters. 



1. Introduction 

The phenomenon of neutrino flavor oscillations, induced by flavor-mass eigenstate mixing, is considered the 
favored solution for solar and atmospheric neutrino data 1 1 1-|5| over a wide range of alternative solutions |6|- 
1141 . These latter mechanisms were considered under the hypothesis that each one of them solely accounts for 
the observed effects. Here we address the possibility of a mixed scenario: one mechanism, the mass-induced 
flavor oscillations, is dominant and a second mechanism is included as sub-dominant: it could be neutrino 
flavor transitions induced by violations of relativity principles, i.e. violation of the Lorentz invariance (VLI) 
or of the equivalence principle (VEP). In this mixed scenario, assuming that neutrinos can be described in 
terms of three distinct bases - flavor, mass and velocity eigenstates - the latter being characterized by different 
maximum attainable velocities (MAVs), and by considering that only two families contribute to the atmospheric 
v violation oscillations, the survival probability can be expressed as 1 101- 1121 

P^u„ = 1- sin 2 26 sin 2 ft (1) 

where the global mixing angle 9 and the term ft are given by: 

28 = arctan(<xi/a2), ft = y (a 2 + a 2 ) . (2) 

The terms a% and a 2 in Eq. (0 contain the relevant physical information 

di = 1.27|Am 2 sin26» m L/ J B + 2-10 18 At)sin26' 1 , J L£;e l ' ) | 

a 2 = 1.27 (Am 2 cos29 m L/E + 2 ■ 10 18 Av cos 26 V LE) , ( } 

where the muon neutrino pathlength L is expressed in km, the neutrino energy E in GeV and the oscillation 
parameters Am 2 = to 2 ™ — m 2 m and At; = iv« — v v v are in eV 2 and c units, respectively. The unconstrained 
phase rj refers to the connection between mass and velocity eigenstates. The whole domain of variability of the 
parameters can be accessed with the requirements Ato 2 > 0, < 9 m < 7r/2, At; > and — 7r/4 < 6 V < 7r/4. 
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The same formalism also applies to violation of the equivalence principle, after substituting Av/2 with the 
adimensional product A7; A7 is the difference of the coupling constants for neutrinos of different types to 
the gravitational potential <j> 1151 . 
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Figure 1. Energy dependence of the i/ M — » v T os- 
cillation probability for mass induced oscillations 
alone (continuous curves), and mass-induced + 
VLI oscillations for Av = 2 x 10~ 25 and dif- 
ferent values of the sin 2 26 v parameter (dashed 
curves). The neutrino pathlength was fixed at 
L = 10 4 km. 
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As shown in II II IT2l 1161 . the most sensitive tests of VLI can be made by analysing the high energy tail 
of atmospheric neutrinos at large pathlengths values. Fig. shows the energy dependence of the — » u T 
oscillation probability as a function of the neutrino energy, for neutrino mass-induced oscillations alone and 
for both mass and VLI-induced oscillations. Notice that for m„ < 1 eV and E v > 100 GeV: 7^ > 10 11 . 



2. Exotic contributions to mass-induced oscillations 

In order to analyse the MACRO data in terms of VLI, we used the subsample of 300 upward-throughgoing 
muons whose energies were estimated via Multiple Coulomb Scattering 1 17 1-| 18 1. Two independent and com- 
plementary analyses have been carried out: one based on the \ 2 criterion and the Feldman and Cousins pre- 
scription 1 19 1, and a second one based on the maximum likelihood technique |20|. 

2.1 x 2 Analysis 

Following the analysis in 1181, we selected a low and a high energy sample by requiring that the reconstructed 
neutrino energy E r v ec should be E r v ec < 30 GeV and E r v ec > 130 GeV. The number of events surviving these 
cuts is Ni ow = 49 and Nhigh = 58; their median energies, estimated via Monte Carlo, are 13 GeV and 204 
GeV (assuming mass-induced oscillations). We then keep the neutrino mass oscillation parameters fixed at 
the values of 1 3 1. The factor e"' is assumed to be real (77 = or it). Then, we scanned the plane of the 2 free 
parameters ( Av, 6 V ) minimizing the \ 2 function comprehensive of statistical and systematic uncertainties [ 20 1 . 
For the Monte Carlo simulation described in 1 1 8 1 the neutrino fluxes in input is given by 1 2 1 1 . The largest rela- 
tive difference of the extreme values of the MC expected ratio Ni ow /Nhigh is 13%. However, in the evaluation 
of the systematic error, the main sources of uncertainties for this ratio (namely the primary cosmic ray spectral 
index and neutrino cross sections) have been separately estimated and their effects added in quadrature (see 1181 
for details): in this work, we use a conservative 16% theoretical systematic error on the ratio Ni ow / Nhi g h-The 
experimental systematic error on the ratio was estimated to be 6%. The inclusion of the VLI effect does not 
improve the \ 2 m an Y point of the (Av, 6 V ) plane, compared to mass-induced oscillations stand-alone, and 
limits on VLI parameters were obtained. The 90% C.L. limits on At; and 6 V , computed with the Feldman 
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and Cousins prescription 1 19 1, are shown by the dashed line in Fig.|2j-left. The energy cuts described above 
were optimized for mass-induced neutrino oscillations. In order to maximize the sensitivity of the analysis for 
VLI induced oscillations, we performed a blind analysis, based only on Monte Carlo events, to determine the 
energy cuts which yield the best performances. The results of this study suggest the cuts E r v ec < 28 GeV and 
E r v ec > 142 GeV; with these cuts the number of events in the real data are N' low = 44 events and N' high = 35 
events. The limits obtained with this selection are shown in Fig.|2l-left by the continuous line. As expected, the 
limits are now more stringent than for the previous choice. In order to understand the dependence of this result 
with respect to the choice of the Am 2 parameter, we varied them around the best-fit point. We found that a 
variation of Am 2 of ±30% moves up/down the upper limit of VLI parameters by at most a factor 2. Finally, 
we computed the limit on Av marginalized with respect to all the other parameters left free to change inside 
the intervals: Am 2 = Am 2 ±30%, 9 m = m ±20%, -7r/4 < 9 V < 7r/4 and any value of the phase 77. We 
obtained the 90% C.L. upper limit \Av\ < 3 x 10~ 25 . 



2.2 Maximum Likelihood Analysis (MLT) 

A different and complementary analysis of VLI contributions to the atmospheric neutrino oscillations was made 
on the MACRO muon data corresponding to parent neutrino energies in the range 25 GeV < E < 75 GeV. 
This energy region is characterized by the best energy reconstruction, and the number of muons satisfying 
this selection is 106. These events are outside the energy ranges used in the analysis discussed in Section 
2.1, and thus the expected sensitivity to VLI (or VEP) contributions to the atmospheric neutrino oscillations 
should be lower; on the other hand, the maximum likelihood technique (MLT) has the advantage to exploit the 
information event-by-event (it is a bin-free approach). Given a specific hypothesis, MLT allows to determine 
the set of parameters a = (Am 2 , 9 m , Av, 9 V )) that maximizes the probability of the realization of the actual 
measurements x = (E, L); it was accomplished by minimizing negative a log-likelihood function C 1221 . We 
have chosen different fixed values of the Am 2 and sin 2 26 m mass-oscillation parameters in |3| and found the 
relative Av and sin 2 29 v that maximize the /(x^; a) function proportional to the probability of realization of a 
given event. Fig. |2]shows the 90% CL upper of the VLI parameter Aw/2 versus the assumed Am 2 values. 
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Figure 2. -left: 90% C.L. upper limits on the Lorentz invariance violation parameter Av versus sin 28 v . Mass induced 
oscillations are assumed in the two-flavor — » v T approximation, with Am 2 = 0.0023 eV 2 and 9 m = 7r/4. The dashed 
line shows the limit obtained with the same selection criteria of Ref. 1 18 1 to define the low and high energy samples; the 
continuous line is the result obtained with the selection criteria optimized for the present analysis (see text), -right: 90% 
CL upper limits on the Av/2 parameter versus Am 2 varying inside the 90% CL MACRO global result [3 ]. 



4 



G. Battistoni et al. 



3. Conclusions 

We have searched for "exotic" contributions to standard mass-induced atmospheric neutrino oscillations arising 
from a possible violation of Lorentz invariance using two different and complementary analyses. The first 
approach uses two sub-sets of events referred to as the low energy and the high energy samples. The mass 
neutrino oscillation parameters have the values of |.3|, and we mapped the evolution of the \ 2 estimator in the 
plane of the VLI parameters, Av and sin 2 26 v . No \ 2 improvement was found, so we applied the Feldman 
Cousins method to determine 90% CL limits on the VLI parameters: \Av\ < 6 x 10~ 24 at sin20„ = and 
| Aw | < 2.5 -f5 x 10~ 26 at sin 26 v = ±1. In terms of the parameter Av alone (marginalization with respect to 
all the other parameters), the VLI parameter bound is (at 90% C.L.) \Av\ < 3 x 10~ 25 . These results may be 
reinterpreted in terms of 90% C.L. limits of parameters connected with violation of the equivalence principle, 
giving the limit |<^A^| < 1.5 x 10~ 25 . The second approach exploits the information contained in a data 
sub-set characterized by intermediate muon energies. It is based on the MLT, and considers the mass neutrino 
oscillation parameters inside the 90% border of the global result |3 1. The obtained 90% CL limit on the At; 
VLI parameter is also around 10 -25 . 
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